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Abstract. We present approaches for the study of fluid-structure interactions subject to thermal 
fluctuations. A mixed mechanical description is utilized combining Eulerian and Lagrangian refer- 
ence frames. We establish general conditions for operators coupling these descriptions. Stochastic 
driving fields for the formalism are derived using principles from statistical mechanics. The stochas- 
tic differential equations of the formalism are found to exhibit significant stiffness in some physical 
regimes. To cope with this issue, we derive reduced stochastic differential equations for several 
physical regimes. We also present stochastic numerical methods for each regime to approximate the 
fluid-structure dynamics and to generate efficiently the required stochastic driving fields. To validate 
the methodology in each regime, we perform analysis of the invariant probability distribution of the 
stochastic dynamics of the fluid-structure formalism. We compare this analysis with results from 
statistical mechanics. To further demonstrate the applicability of the methodology, we perform com- 
putational studies for spherical particles having translational and rotational degrees of freedom. We 
compare these studies with results from fluid mechanics. The presented approach provides for fluid- 
structure systems a set of rather general computational methods for treating consistently structure 
mechanics, hydrodynamic coupling, and thermal fluctuations. 
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1. Introduction. The development of analytic and computational approaches 
for the study of fluid-structure interactions has a rich history. Motivations for past 
work in this area include the study of aerodynamic oscillations induced in airplane 
wings and propellers [22, 24], the study of animal locomotion including swimming and 
insect flight [44, 38, 51], and the study of physiological problems such as blood flow 
through heart valves [30, 49, 27]. A central challenge in work on these applications 
has been to develop descriptions which capture essential features of the fluid struc- 
ture interactions while introducing approximations which facilitate analysis and the 
development of tractable numerical methods [22, 45]. Many such challenges remain 
and this area of research is still very active [30, 13, 45, 12]. Recent scientific and 
technological advances motivate the study of fluid-structure interactions in new phys- 
ical regimes often involving very small length scales [57, 64, 15, 46]. At sufficiently 
small length scales thermal fluctuations play an important role and pose additional 
challenges in the study of fluid-structure systems. 

Significant past work has been done on the formulation of descriptions of fluid- 
structure interactions subject to thermal fluctuations. Many of these analytic and nu- 
merical approaches originate from the polymer physics community [19, 23, 52, 11]. To 
obtain descriptions tractable for analysis and numerical simulation, these approaches 
typically place an emphasis on approximations which retain only the structure degrees 
of freedom. This often results in significant simplifications in the descriptions and in 
significant computational savings. This eliminates the many degrees of freedom as- 
sociated with the fluid and avoids having to resolve the potentially intricate and stiff 
stochastic dynamics of the fluid. These approaches have worked especially well for 
the study of bulk phenomena in free solution and the study of complex fluids and soft 
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materials [52, 19, 37]. 

Recent applications arising in the sciences and in technological fields present situ- 
ations in which resolving the dynamics of the fluid may be important and even advan- 
tageous both for modeling and computation. This includes modeling the spectroscopic 
responses of biological materials [65, 28, 43], studying transport in microfluidic and 
nanofmidic devices [57, 47], and investigating dynamics in biological systems [2, 17]. 
There are also other motivations for representing the fluid explicitly and resolving its 
stochastic dynamics. This includes the development of hybrid fluid-particle models in 
which thermal fluctuations mediate important effects when coupling continuum and 
particle descriptions [18, 20], the study of hydrodynamic coupling and diffusion in the 
vicinity of surfaces having complicated geometries [57], and the study of systems in 
which there are many interacting mechanical structures [7, 50, 49]. To facilitate the 
development of methods for studying such phenomena in fluid-structure systems, we 
present a rather general formalism which captures essential features of the coupled 
stochastic dynamics of the fluid and structures. 

To model the fluid-structure system, a mechanical description is utilized involving 
both Eulerian and Lagrangian reference frames. Such mixed descriptions arise rather 
naturally, since it is often convenient to describe the structure configurations in a 
Lagrangian reference frame while it is convenient to describe the fluid in an Eulerian 
reference frame. In practice, this presents a number of challenges for analysis and 
numerical studies. A central issue concerns how to couple the descriptions to represent 
accurately the fluid-structure interactions, while obtaining a coupled description which 
can be treated efficiently by numerical methods. Another important issue concerns 
how to account properly for thermal fluctuations in such approximate descriptions. 
This must be done carefully to be consistent with statistical mechanics. A third 
issue concerns the development of efficient computational methods. This requires 
discretizations of the stochastic differential equations and the development of efficient 
methods for numerical integration and stochastic field generation. 

We present a set of approaches to address these issues. The formalism and general 
conditions for the operators which couple the Eulerian and Lagrangian descriptions 
arc presented in Section 2. We discuss simplified descriptions of the fluid-structure 
system in different physical regimes in Section 3. A derivation of the stochastic driv- 
ing fields used to represent the thermal fluctuations is also presented in Section 3. 
Stochastic numerical methods are discussed for the approximation of the stochastic 
dynamics and generation of stochastic fields in Sections 4. To validate the method- 
ology, we perform in each regime analysis of the invariant probability distribution of 
the stochastic dynamics of the fluid-structure formalism. We compare this analysis 
with results from statistical mechanics in Section 5. To demonstrate the applicability 
of the methodology, we perform computational studies for spherical particles having 
translational and rotational degrees of freedom. We compare these computational 
studies with results from fluid mechanics in Section 6. 

It should be mentioned that related computational methods have been intro- 
duced for the study of fluid-structure interactions [50, 3, 10, 12, 63, 45, 67, 33, 41]. 
In recent papers, significant work also has been done toward incorporating the role of 
thermal fluctuations [6, 14, 9, 21, 7]. This includes the Stochastic Immersed Boundary 
Method [6], Fluctuating Immersed Material Dynamics [14], Computational Fluctuat- 
ing Fluid Dynamics [9, 21], and Accelerated Stokesian Dynamics [7]. The formalism 
presented here can be regarded in part as a generalization of these approaches. It is 
expected that many of the presented results can be applied to further justify and val- 
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Fig. 2.1. The description of the fluid- structure system utilizes both Eulerian and Lagrangian 
reference frames. The structure mechanics are often most naturally described using a Lagrangian 
reference frame. The fluid mechanics are often most naturally described using an Eulerian reference 
frame. The mapping X(q) relates the Lagrangian reference frame to the Eulerian reference frame. 
The operator T prescribes how structures are to be coupled to the fluid. The operator A prescribes 
how the fluid is to be coupled to the structures. A variety of fluid-structure interactions can be 
represented in this way. This includes rigid and deformable bodies, membrane structures, polymeric 
structures, or point particles. 



idate these methods and to provide further extensions. The formalism presented here 
provides a rather general framework for the development of computational methods for 
applications requiring a consistent treatment of structure mechanics, hydrodynamic 
coupling, and thermal fluctuations. 

2. Summary of the Stochastic Eulerian Lagrangian Method . We sum- 
marize here the Stochastic Eulerian Lagrangian Method, abbreviated as SELM. We 
present the general formalism and a number of alternative descriptions of the fluid- 
structure system. In many situations the stochastic differential equations for the full 
fluid-structure dynamics exhibits stiffness. To cope with this issue and to develop 
efficient numerical methods, simplified descriptions are discussed for various physi- 
cal regimes. A more detailed discussion and derivation of SELM and the reduced 
equations in each of the physical regimes is given in Section 3. 

To study the dynamics of fluid-structure interactions in the presence of thermal 
fluctuations, we utilize a mechanical description involving Eulerian and Lagrangian 
reference frames. Such mixed descriptions arise rather naturally, since it is often 
convenient to describe the structure configurations in a Lagrangian reference frame 
while it is convenient to describe the fluid in an Eulerian reference frame. In princi- 
ple more general descriptions using other reference frames could also be considered. 
Descriptions for fluid-structure systems having these features can be described rather 
generally by the following dynamic equations 

(2.1) =£u + A[T(v-ru)] + A + f thm 
at 

(2.2) = -T (v - Tu) - V X $[X] + C + F thm 
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The u denotes the velocity of the fluid, p the uniform fluid density. The X denotes 
the configuration of the structure and v the velocity of the structure. The mass of 
the structure is denoted by m. To simplify the presentation we treat here only the 
case when p and m are constant, but with some modifications these could also be 
treated as variable. The A, Q are Lagrange multipliers for imposed constraints, such 
as incompressibility of the fluid or a rigid body constraint of a structure. The operator 
C is used to account for dissipation in the fluid, such as associated with Newtonian 
fluid stresses [1]. To account for how the fluid and structures are coupled, a few 
general operators are introduced, r,T,A. 

The linear operators T, A, T are used to model the fluid-structure coupling. The 
F operator describes how a structure depends on the fluid flow while — T is a negative 
definite dissipative operator describing the viscous interactions coupling the structure 
to the fluid. We assume throughout that this dissipative operator is symmetric, 
T = T T . The linear operator A is used to attribute a spatial location for the viscous 
interactions between the structure and fluid. The linear operators are assumed to 
have dependence only on the configuration degrees of freedom V = T[X], A = A[X]. 
We assume further that T does not have any dependence on X. 

To account for the mechanics of structures, 3>[X] denotes the potential energy of 
the configuration X. The total energy associated with this fluid-structure system is 
given by 



The first two terms give the kinetic energy of the fluid and structures. The last term 
gives the potential energy of the structures. 

As we shall discuss, it is natural to consider coupling operators A and T which 
are adjoint in the sense 



for any u and v. The S and Q denote the spaces used to parameterize respectively the 
structures and the fluid. We denote such an adjoint by A = 1^ or V = A^. This adjoint 
condition can be shown to have the important consequence that the fluid-structure 
coupling conserves energy when T — > oo in the inviscid and zero temperature limit. 

In practice, the conditions discussed above can be relaxed somewhat. For our 
present purposes these conditions help simplify the presentation. Each of these oper- 
ators will be discussed in more detail below. 

To account for thermal fluctuations, a random force density f thm is introduced 
in the fluid equations and F thm in the structure equations. These account for spon- 
taneous changes in the system momentum which occurs as a result of the influence 
of unresolved microscopic degrees of freedom and unresolved events occurring in the 
fluid and in the fluid-structure interactions. 

The thermal fluctuations consistent with the form of the total energy and relax- 
ation dynamics of the system are taken into account by the introduction of stochastic 
driving fields in the momentum equations of the fluid and structures. The stochastic 
driving fields are taken to be Gaussian processes with mean zero and with ^-correlation 
in time [54] . By the fluctuation-dissipation principle [54] these have covariances given 



(2.4) 




(2.5) 




by 



(2.6) 



{2k B T) (C - ATT) S(t - s) 
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(2.7) 
(2.8) 



(F th Us)FlJt)) = (2k B T)TS(t-s) 
(f thm ( s )Fl,(<)) = - (2k B T) AT S(t - s). 



We have used that T = and T = T T . We remark that the notation gh T which is 
used for the covariance operators should be interpreted as the tensor product. This 
notation is meant to suggest the analogue to the outer-product operation which holds 
in the discrete setting [5] . A more detailed discussion and derivation of the thermal 
fluctuations is given in Section 3. 

It is important to mention that some care must be taken when using the above 
formalism in practice and when choosing operators. An important issue concerns 
the treatment of the material derivative of the fluid, du/dt = du/dt + u • Vu. For 
stochastic systems the field u is often highly irregular and not defined in a point-wise 
sense, but rather only in the sense of a generalized function (distribution) [16, 40] . This 
presents issues in how to define the non-linear term arising in the material derivative, 
which appears to require point-wise values of u. For such irregular velocity fields, 
this also calls into question the applicability of the theorems typically used to derive 
the differential equations from the conservation laws. For instance, for such velocity 
fields the fluid material body may no longer exhibit smooth deformations over time. 

There are a number of ways to deal with this issue. The first is to consider 
a regularization of the fluid stresses, which are typically the source of irregularity, 
see equation 2.6. This can be motivated by the fact that the fluid stress tensors 
typically considered in continuum mechanics are expected to become inaccurate at 
molecular length-scales. Ideally, from molecular models of the fluid the small-length 
scale (large wave-number) responses of the fluid could be determined and provide 
a justified regularization. For instance, this could provide an alternative to using 
responses based on Newtonian stresses for all length-scales. For the SELM formalism, 
this would simply correspond to using for L an alternative to the dissipative operator 
based on Newtonian stresses. The second more easily implemented approach is simply 
to work with the linearized material derivative, which still retains many of the essential 
features of the fluid dynamics and is useful for many applications [6]. 

In this initial presentation of SELM, we shall take the latter approach and treat 
du/dt = du/dt. This provides a rather general description of fluid-structure sys- 
tems which incorporate the role of thermal fluctuations. From this initial formalism 
of SELM, we shall derive a number of simplified descriptions for various physical 
regimes. These simplified descriptions for each regime tend to yield less stiff differ- 
ential equations and have other features making them useful in the development of 
efficient stochastic numerical methods for the formalism. 

2.1. Regime I. We now consider the regime in which the full dynamics of the 
fluid-structure system are retained but reformulated in terms of a field describing the 
total momentum of the fluid-structure system at a given spatial location. This de- 
scription is more convenient to work with in practice since it results in simplifications 
in the stochastic driving fields. For this purpose we define 



The operator A is used to give the distribution in space of the momentum associated 
with the structures for given configuration X(i). Using this approach, the fluid- 
structure dynamics are described by 



(2.9) 



p(x,i) = pu{x,t) + A[mv(t)](x). 



(2.10) 



dp 

It 



Cu + A[-V X $(X)] + (V x A[mv]) ■ v + A + g, 



.th] 
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(2.11) m^ = -T(v-ru)-V x $(X) + C + F t 



(2.12) 



dt 
~dt 



where u = p -1 (p — A[rav]) and g thm = f thm + A[F thm ]. The third term in the first 
equation arises from the dependence of A on the configuration of the structures, 
A[mv] = (A[X])[rav]. The Lagrange multipliers for imposed constraints are denoted 
by A, C- For the constraints, we use rather liberally the notation with the Lagrange 
multipliers denoted here not necessarily assumed to be equal to the previous definition. 
The stochastic driving fields are again Gaussian with mean zero and ^-correlation in 
time [54] . The stochastic driving fields have the covariance structure given by 

(2-13) (Bt.un(«)gL(*)) = - ( 2fc sT) C S(t - s) 

(2-14) (F thm (s)FL W) - (2k B T) T S(t - s) 

(2-15) (g thm ( S )FLW) = 0. 

This formulation has the convenient feature that the stochastic driving fields become 
independent. This is a consequence of using the field for the total momentum for 
which the dissipative exchange of momentum between the fluid and structure no 
longer arises. In the equations for the total momentum, the only source of dissipa- 
tion remaining occurs from the stresses of the fluid. This approach simplifies the 
effort required to generate numerically the stochastic driving fields and will be used 
throughout. 

2.2. Regime II. We now consider a regime in which the formalism can be sim- 
plified significantly. In many situations, inertial effects often play a relatively minor 
role in the structure dynamics as a consequence of the small mass of the structure rela- 
tive to the displaced fluid or as a consequence of viscosity of the solvent fluid [31, 42]. 
In such a regime, the relatively rapid dynamics associated with the momentum of 
the structures often presents a source of stiffness in numerical calculations. To cope 
with these issues we consider a reduction of the stochastic dynamics of the system in 
which the structure momentum is eliminated from the description. In particular, we 
consider the regime in which m -C p( 3 ■ The i denotes a length-scale characteristic of 
the size of the immersed structure and associated flow field of the fluid. In the limit 
m — > 0, the fluid-structure dynamics are governed by 

(2.16) ^ = p-'Cp + A[-V X $(X)] + (V x ■ A) k B T + A + g thm 

(2.17) ^ - p-'Tp + T- 1 [-V x <f(X)] + C + G thm . 

In the notation Vx ■ A = Tr[VxA]. This term arises from the thermal fluctuations 
associated with the momentum of the structures, which have been eliminated from 
the description. This term plays an important role in the system when the phase- 
space dynamics of (p, X) has an associated vector field which is compressible. When 
considering the Liouville equation on phase-space, this term accounts for local changes 
of the phase-space volume which occurs as the configuration of the structure changes 
under the dynamics of the reduced description [62]. For a more detailed discussion 
see Section 5.2. The stochastic driving fields g th m,G thm are again Gaussian with 
mean zero and with ^-correlation in time [54]. The stochastic driving fields have the 
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covariance structure given by 

(2-18) (g thm (s)g t T hm («)) = - (2k B T) C S(t - s) 

(2.19) (G ttm ( S )Gl(t)) = (2k B T) IT 1 S(t - a) 

(2-20) (gtaWGLW) = 0. 

A more detailed discussion and derivation of the equations in this regime is given in 
Section 3.2. 

2.3. Regime III. The description of the fluid-structure system can be further 
simplified by considering the viscous coupling between the fluid and structures in the 
limit T — > oo. In this case, the fluid-structure dynamics are given by 

(2.21) ^ = p^Cp + A[-V X *(X)] + (Vx • A) k B T + X + g thm 
at 

(2.22) ^ = p^Tp 

(2.23) ( gthm (s)gL(0> = - (2k B T)C6(t-s). 

In the notation Vx • A = Tr[VxA]. A more detailed discussion and derivation of the 
equations in this regime in given in Section 3.3. 

2.4. Regime IV. The description of the fluid-structure system can be further 
simplified by considering for the fluid the viscous limit in which fj, — > oo. In this 
regime only the structure dynamics remain and can be shown to be given by 

(2.24) ^ = # SELM [-V X $(X)] + (V x • H SEhM )k B T + h thm 

(2.25) ff SELM = T(-p£)- 1 A 

(2.26) (h thm (s)h t T hm (t)) = (2k B T) ff SELM S(t - s). 

The p denotes a projection operator imposing constraints, such as incompressibility. 
The adjoint property A = and symmetry of pC yields an operator i? SBLM which is 
symmetric. A more detailed discussion and derivation of the equations in this regime 
is given in Section 3.4. 

2.5. Summary. This gives an overview of the SELM formalism and the asso- 
ciated stochastic differential equations. We remark that each of these regimes were 
motivated by a rather specific limit. Non-dimensional analysis of the equations can 
also be carried out and other limits considered to motivate working with such reduced 
equations. We discuss in more detail the derivation of the reduced equations in each 
regime in Section 3. We discuss how specific stochastic numerical methods can be 
developed for the SELM formalism in Section 4. We discuss applications and how the 
SELM formalism can be used in practice in Section 6. 

3. Derivations for the Stochastic Eulerian Lagrangian Method. We now 

discuss formal derivations to motivate the stochastic differential equations used in 
each of the physical regimes. For this purpose, we do not present the most general 
derivation of the equations. For brevity, we make simplifying assumptions when 
convenient. 

In the initial formulation of SELM, the fluid-structure system is described by 

(3.1) /0 ^ = £u + A[T(v-ru)] + A + f thin 

at 
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(3.2) 
(3.3) 



dv 

'T'—r 
dt 

dX 

~dt 



T (v - ru) - V X $(X) + C + F thn 



The notation and operators appearing in these equations has been discussed in de- 
tail in Section 2. For these equations, we focus primarily on the motivation for the 
stochastic driving fields used for the fluid-structure system. 

For the thermal fluctuations of the system, we assume Gaussian random fields 
with mean zero and 5-correlated in time. For such stochastic fields, the central chal- 
lenge is to determine an appropriate covariance structure. For this purpose, we use the 
fluctuation-dissipation principle of statistical mechanics [54, 36]. For linear stochastic 
differential equations of the form 



(3.4) 



dZ t = LZ t dt + QdB t 



the fluctuation-dissipation principle can be expressed as 
(3.5) G = QQ T = —(LC) — (LC) T . 

This relates the equilibrium covariance structure C of the system to the covariance 
structure G of the stochastic driving field. The operator L accounts for the dissipative 
dynamics of the system. For the equations 3.1 - 3.3, the dissipative operators only 
appear in the momentum equations. This can be shown to have the consequence that 
there is no thermal forcing in the equation for X(i), this will also be confirmed in 
Section 5.1. To simplify the presentation, we do not represent explicitly the stochastic 
dynamics of the structure configuration X. 

For the fluid-structure system it is convenient to work with the stochastic driving 
fields by defining 



(3.6) 



q = [p ^tto 



The field q formally is given by q = QdB t /dt and determined by the covariance 
structure G = QQ T . This covariance structure is determined by the fluctuation- 
dissipation principle expressed in equation 3.5 with 



(3.7) 
(3.8) 



L 



C 



p- 1 (C - ATT) 

m^rr 





P 



L AT 



Tl 



o 

m ~ 



l k R TI 



The I denotes the identity operator. The covariance C was obtained by considering 
the fluctuations at equilibrium. The covariance C is easily found since the Gibbs- 
Boltzmann distribution is a Gaussian with formal density ^(u, v) = ~ exp [— E/kgT]. 
The Zq is the normalization constant for 'J'. The energy is given by equation 2.4. For 
this purpose, we need only consider the energy E in the case when $ = 0. This gives 
the covariance structure 



(3.9) 



G = (2k B T) 



-p~ 2 (C — ATr) -m-ip- 1 AT 



To obtain this result we use that r = At and T = T"i". From the definition of 
q, it is found the covariance of the stochastic driving fields of SELM are given by 
equations 2.6- 2.8. This provides a description of the thermal fluctuations in the 
fluid-structure system. 
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3.1. Regime I. It is convenient to reformulate the description of the fluid- 
structure system in terms of a field for the total momentum of the system associated 
with spatial location x. For this purpose we define 

(3.10) P(x,£) = pu(x,t) + A[mv(i)](x). 

The operator A is used to give the distribution in space of the momentum associated 
with the structures. Using this approach, the fluid-structure dynamics are described 

by 

(3.11) d ^ t =Cn + A h v x$(X)] + (V X AH) ■ v + A + gthm 

(3.12) m^r = -T (v - Tu) - V X $(X) + ( + F thm 



(3.13) 



dt 
rfX 

~dt 



where u = p 1 (p — A[mv]) and g thm = f thm + A[F thm ]. The third term in the first 
equation arises from the dependence of A on the configuration of the structures, 
A[mv(i)] = (A[X])[mv(t)]. 

The thermal fluctuations are taken into account by two stochastic fields g thm and 
F thm . The covariance of g thm is obtained from 

(3.14) (g thm gL) = (fiX) + (f t .»FL A T ) + (AF thm f t L) + (AF thm Fl A T ) 
= (2k B T) (-£ + ATT - ATA T - ATA T + ATA T ) 
= -(2k B T)C. 

This makes use of the adjoint property of the coupling operators A^ = T. 

One particularly convenient feature of this reformulation is that the stochastic 
driving field F thm and g thm become independent. This can be seen as follows 

(3-15) (g thm F t T hm ) = (f thm FL) + (AF thm Ff hm ) 

= (2fc B T)(-AT + AT) = 0. 

This decoupling of the stochastic driving fields greatly reduces the computational ef- 
fort to generate the fields with the required covariance structure. This shows the co- 
variance structure of the stochastic driving fields of SELM are given by equations 2.13— 
2.15. 

3.2. Regime II. In many situations, inertial effects often play a relatively minor 
role in the structure dynamics as a consequence of the small mass of the structure 
relative to the displaced fluid or as a consequence of the large viscosity of the solvent 
fluid [31, 42]. We consider the regime in which m -C p£ 3 , as discussed in Section 2.2. 
We shall derive formally reduced stochastic equations in the limit m — > 0. 

For this purpose, we focus primarily on the dynamics of the velocity of the struc- 
tures v(<). This can be expressed using the notation of Ito Stochastic Differential 
Equations [48] as 

(3.16) dV t = -mT l T (V t - Tu + T^Vx^KX)) dt + m- 1 (2k B TT) 1/2 dB t . 

The B t denotes throughout the standard Brownian motion on M. N [48]. To simplify 
the presentation, we consider the case when £ = 0. We expect similar results to 
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hold more generally. Treating the other degrees of freedom as fixed, we can solve 
equation 3.16 using Ito's Lemma [48]. The stationary behavior of this stochastic 
process can be expressed as 



/t 
e -(t- s )m 1 v i (2jt B n) 1 / 2 dB s 
-oo 

(3.18) fi = ru-T- 1 V x $(X). 

As a result of the integrand being deterministic in the Ito Integral, the V t is a Gaussian 
process. This has the consequence that the statistics of the process V t are completely 
determined by its mean and covariance functions. The mean of the process is given 
at each time by 

(3.19) n(t) = (V t ) = fi Q . 

The covariance function can be computed using the Ito Isometry [48] to obtain 

(3.20) 0(|r|) = ((V t+T - n )(V t - Mo ) T ) = kBTm^e-^^l. 
In the limit m — >■ this can be expressed as 

(3.21) 4>{\t - s\) = 2k B TT- 



2 



2k B TT- 1 6(t- s). 



We have used formally ^Ae A ' r ' — > 5(t) as A — > oo. This suggests the following 
approximation for the velocity of the structures in equations 2.11 and 2.12. 

(3.22) v(t) -> Tu - T- 1 V X *(X) + (2k B TT- 1 f /2 

To approximate the term (VxA[mv]) • v in the limit m — > appearing in equa- 
tion 2.10, a different approach is required. For this purpose, we consider the process 
R t — /(Vt) = mVtVf- By Ito's Lemma this satisfies the stochastic differential 
equation 

(3.23) dR t = V/(V t )dVt + \d\J V 2 /(V t )dV t 

with the formal substitutions dtdRt = = dtdt, dB t dSj = Xdt. To simplify the 
discussion we consider the case when T = jX. We expect similar results can be 
obtained more generally. In this case 

(3.24) dR t = -2m- 1 7 (^ - k B TI - ^(mV t ^ + fi (mV t ) T ))dt 

+ m- 1 {2k B T 1 ) 1/2 (mVtdBf + dB t (mV t ) T ) . 

From the form of the energy given by equation 2.4, the structure momentum mV f 
has equilibrium distribution ^(mV) = exp [— (rriV) 2 /2mk B T] . This gives a Gaussian 
with mean and variance 

(3.25) (mV t ) = 

(3.26) ((mV t )(mV() T ) = mk B Tl. 
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By reasoning similar to the arguments used to approximate V t , this suggests the 
terms involving mV t do not make a contribution in the m — > limit. This suggests 
to leading order we have 

(3.27) dR t = -2m~ l -f(R t - k B Tl)dt. 
This suggests the approximation in the limit m — > 

(3.28) R t = mV t V t T -> k B T2. 
By substituting this result for v in equation 2.10, we have 

(3.29) (V X AH) ■ v -> Tr[V x A](fc B T) = (V x • A)(A B T). 

This establishes formally the reduced SELM description when the mass of the 
structures becomes negligible. It should be mentioned these approximations can be 
established more rigorously using a perturbation analysis of the Kolomogorov Equa- 
tions associated with the stochastic processes [35, 25]. It should also be mentioned 
that other limits can be considered in which additional drift and stochastic terms arise 
in the reduced momentum equations. This will be the focus of another paper. 

3.3. Regime III. The SELM stochastic equations can be further reduced if the 
effective viscous interactions between the structure and fluid are assumed to become 
very strong. This corresponds to approximating the reduced stochastic equations of 
Regime II in the formal limit T — > oo. By this notation, we mean that all eigenvalues 
of the symmetric operator T uniformly tend to infinity. In this formal limit the terms 
involving T _1 are expected to no longer make a contribution to the dynamics. This 
motivates the reduced stochastic equations 2.21- 2.23. 

3.4. Regime IV. The description of the fluid-structure system can be further 
simplified by considering for the fluid the viscous limit in which \i — > oo . In this regime 
the fluid adopts a quasi-steady-state behavior with respect to the configuration of the 
structures and the forces acting on the fluid. In this limit only the structure dynamics 
remain. By approximating u(i) using arguments similar to those used in Regime II 
for approximating V t , we can derive the reduced stochastic equations 2.24- 2.26. 

4. Computational Methodology. We now discuss briefly numerical methods 
for the SELM formalism. For concreteness we consider the specific case in which the 
fluid is Newtonian and incompressible. For now, the other operators of the SELM 
formalism will be treated rather generally. This case corresponds to the dissipative 
operator for the fluid 

(4.1) Cu = /zAu. 

The A denotes the Laplacian Au = d xx u + d yy u + d zz u. The incompressibility of the 
fluid corresponds to the constraint 

(4.2) V • u = 0. 

This is imposed by the Lagrange multiplier A. By the Hodge Decomposition, A is 
given by the gradient of a function p with A = — Vp. The p can be interpreted as the 
local pressure of the fluid. 

A variety of methods could be used in practice to discretize the SELM formal- 
ism, such as Finite Difference Methods, Spectral Methods, and Finite Element Meth- 
ods [29, 60, 59]. We present here discretizations based on Finite Difference Methods. 
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4.1. Numerical Semi-Discretizations for Incompressible Newtonian Fluid. 

The Laplacian will be approximated by central differences on a uniform periodic lat- 
tice by 



(4-3) [£u] m = £ 



j_„„. - 2u, 



Ax 2 

3 = 1 

The m = (mi, m2, m^) denotes the index of the lattice site. The e^ denotes the 
standard basis vector in three dimensions. The incompressibility of the fluid will be 
approximated by imposing the constraint 

(4-4) [D ■ u] m = ± ■■ 2A ; 4 - ■ 

j=i 

The superscripts denote the vector component. In practice, this will be imposed by 
computing the projection of a vector u* to the sub-space {u £ M. 3N \ D ■ u = 0}, where 
N is the total number of lattice sites. We denote this projection operation by 

(4.5) u=pu*. 

The semi-discretized equations for SELM to be used in practice are 

(4.6) ^ = Lu + A[-V x *] + (V x A[mv]) • v + A + g thm 

(4.7) ^ = -T[v-ru]+F thm 
<iX 

The component u m = p~ 1 (p m — A[mv] m ). Each of the operators now appearing 
are understood to be discretized. We discuss specific discretizations for T and A in 
Section 6. To obtain the Lagrange multiplier A which imposes incompressibility we 
use the projection operator and 

(4.9) A = -{X - p) (Lu + T[v - Tu] + f thm ) 

In this expression, we let f thm = g thm — A[F thm ] for the particular realized values of 
the fields g thm and F thm . 

We remark that in fact the semi-discretized equations of the SELM formalism 
in this regime can also be given in terms of u directly, which may provide a simpler 
approach in practice. The identity f thm = g thm — A[F thm ] could be used to efficiently 
generate the required stochastic driving fields in the equations for u. We present the 
reformulation here, since it more directly suggests the semi-discretized equations to 
be used for the reduced stochastic equations. 

For this semi-discretization, we consider a total energy for the system given by 

(4.10) E[u, v, X] = |£ |u(x m )| 2 AxL + ||v| 2 + 4>[X]. 



This is useful in formulating an adjoint condition 2.5 for the semi-discretized system. 
This can be derived by considering the requirements on the coupling operators T 
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and A which ensure the energy is conserved when T — > oo in the inviscid and zero 
temperature limit. 

To obtain appropriate behaviors for the thermal fluctuations, it is important to 
develop stochastic driving fields which are tailored to the specific semi-discretizations 
used in the numerical methods. Once the stochastic driving fields are determined, 
which is the subject of the next section, the equations can be integrated in time using 
traditional methods for SDEs, such as the Euler-Maruyama Method or a Stochastic 
Runge-Kutta Method [34]. More sophisticated integrators in time can also be devel- 
oped to cope with sources of stiffness, but are beyond the scope of this paper [6]. For 
each of the reduced equations, similar semi-discretizations can be developed as the 
one presented above. 

4.2. Stochastic Driving Fields for Semi-Discretizations. To obtain behav- 
iors consistent with statistical mechanics, it is important stochastic driving fields be 
used which are tailored to the specific numerical discretization employed [6, 21, 5]. 
To ensure consistency with statistical mechanics, we will again use the fluctuation- 
dissipation principle but now apply it to the semi-discretized equations. For each 
regime, we then discuss the important issues arising in practice concerning the effi- 
cient generation of these stochastic driving fields. 

4.3. Regime I. To obtain the covariance structure for this regime, we apply the 
fluctuation-dissipation principle as expressed in equation 3.5 to the semi-discretized 
equations 4.6- 4.8. This gives the covariance 



(4.11) G = -2LC = (2k B T) 



-p- 2 Ax- 3 L 
m- 2 T 





The factor of Ax~ 3 arises from the form of the energy for the discretized system which 
gives covariance for the equilibrium fluctuations of the total momentum p^ 1 Ax~ 3 ksT, 
see equation 4.10. In practice, achieving the covariance associated with the dissipative 
operator of the fluid L is typically the most challenging to generate efficiently. This 
arises from the large number N of lattice sites in the discretization. 

One approach is to determine a factor Q such that the block G pp = QQ T , sub- 
scripts indicate block entry of the matrix. The required random field with covariance 
Gp.p is then given by g = Q£, where £ is the uncorrelated Gaussian field with the 
covariance structure X. For the discretization used on the uniform periodic mesh, 
the matrices L and C are cyclic [58]. This has the important consequence that they 
are both diagonalizable in the discrete Fourier basis of the lattice. As a result, the 
field f thm can be generated using the Fast Fourier Transform (FFT) with at most 
0(N \og(N)) computational steps. In fact, in this special case of the discretization, 
"random fluxes" at the cell faces can be used to generate the field in O(N) compu- 
tational steps [5]. Other approaches can be used to generate the random fields on 
non-periodic meshes and on multi- level meshes, see [4, 5]. 

4.4. Regime II. The covariance structure can be found using an approach sim- 
ilar to the one presented in Section 4.3. This gives 



(4.12) G = {2k B T) 



-p- z Ax- A L 



By factoring the covariance matrix in the Fourier basis, the field can be generated 
using FFTs in at most 0(N\og(N)) computational steps. 
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4.5. Regime III. The covariance structure can be found using an approach 
similar to the one presented in Section 4.3. This gives 



(4.13) G = (2k B T) 



-p- 2 Ax~ 3 L 




By factoring the covariance matrix in the Fourier basis, the field can be generated 
using FFT in at most 0(N \og(N)) computational steps. 

4.6. Regime IV. This regime differs from the others since the fluid momentum 
and structure momentum are both no longer represented explicitly. Spontaneous 
changes in the momentum of the system were the primary source of fluctuations 
in the configuration of the structures in the other regimes. While the momentum 
is no longer represented explicitly, we can none-the-less use a discretization of the 
momentum equations to generate efficiently the random fields required in the over- 
damped dynamics. This is done by expressing the covariance of the stochastic driving 
field as 

(4.14) G = (2fc B T)ff SELM = (2k B T) (Tp(-L)' 1 p T T T ) . 

This makes use of A = T T and properties of the specific discretized operators L and 
p. In particular, commutativity pL = Lp and the projection operator properties 
p 2 = p, p — p T . Let U be a factor so that UU T = —L^ 1 . Using this factor we can 
express the covariance as 

(4.15) G = U/2k^fTpU) (j2k B TYpU 

From this expression a matrix square-root of G is readily obtained, Q = ^/2kBTTpU . 

We remark this is different than the Cholesky factor obtained for G which is 
required to be lower triangular [61, 58]. Obtaining such a factor by Cholesky factor- 
ization would cost (3(M 3 ), where M is the number of structure degrees of freedom. 
For the current discretization considered, the operators L and p are diagonalizable in 
Fourier space. This has the consequence that the action of the operators U and p can 
be computed using FFTs with a cost of 0(N log(iV)). The N is the number of lattice 
sites used to discretize L. The stochastic driving field is computed from h = Q£. This 
allows for the stochastic driving field to be generated in 0(N\og(N) + M) computa- 
tional steps, assuming the action A can be compute in 0{M) steps. This is in contrast 
to using the often non-sparse matrix arising from Cholesky factorization which gener- 
ates the stochastic field with a cost of 0(M 2 ). We remark that this approach shares 
some similarities with the method proposed in [7, 56]. Other methods based on split- 
tings or multigrid can also be utilized to efficiently generate stochastic fields with this 
required covariance structure, see [4, 5]. 

5. Equilibrium Statistical Mechanics of SELM Dynamics. We now dis- 
cuss how the SELM formalism and the presented numerical methods capture the equi- 
librium statistical mechanics of the fluid-structure system. This is done through an 
analysis of the invariant probability distribution of the stochastic dynamics. For the 
fluid-structure systems considered, the appropriate probability distribution is given 
by the Gibbs-Boltzmann distribution 

(5.1) $ OB ( Z ) = iexp[-£( z )/fc B T]. 
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The z is the state of the system, E is the energy, ks is Boltzmann's constant, T is the 
system temperature, and Z is a normalization constant for the distribution [54]. We 
show this Gibbs-Boltzmann distribution is the equilibrium distribution of both the 
full stochastic dynamics and the reduced stochastic dynamics in each physical regime. 

We present here both a verification of the invariance of the Gibbs-Boltzmann 
distribution for the general formalism and for numerical discretizations of the formal- 
ism. The verification is rather formal for the undiscretized formalism given technical 
issues which would need to be addressed for such an infinite dimensional dynamical 
system. However, the verification is rigorous for the semi-discretization of the for- 
malism, which yields a finite dimensional dynamical system. The latter is likely the 
most relevant case in practice. Given the nearly identical calculations involved in the 
verification for the general formalism and its semi-discretizations, we use a notation 
in which the key differences between the two cases primarily arise in the definition 
of the energy. In particular, the energy is understood to be given by equation 2.4 
when considering the general SELM formalism and equation 4.10 when considering 
semi-discretizations. 

5.1. Regime I. The stochastic dynamics given by equations 2.10- 2.12 is a 
change-of-variable of the full stochastic dynamics of the SELM formalism given by 
equations 2.1- 2.3. Thus verifying the invariance using the reformulated description 
is also applicable to equations 2.1- 2.3 and vice versa. To verify the invariance in 
the other regimes, it is convenient to work with the reformulated description given 
for Regime I. The energy associated with the reformulated description is given by 

(5.2) £[p, v, n = Y p j n IP(y) - A[mv](y)| 2 dy + ||v| 2 + $[X]. 
The energy associated with the semi-discretization is 

(5.3) E[p, v,X] = l^ |p(x m ) - A[mv] m | 2 Ax 3 m + ™ |v| 2 + $[X]. 

m 

The probability density ^(p, v, X, t) for the current state of the system under the 
SELM dynamics is governed by the Fokker-Planck equation 

(5.4) -— = —V ■ J 



dt 



with probability flux 
(5.5) J = 



C + A + V X A • v + A 

-T- V x $ + C 
v 



The covariance operator G is associated with the Gaussian field g = [g thm , F thm , 0] T 
by (g( s )g T (^)) = G6(t—s). In this regime, G is given by equation 2.13 or 4.11. In the 
notation [V • G(z)]j = d Zj Gij(z) with the summation convention for repeated indices. 
To simplify the notation we have suppressed denoting the specific functions on which 
each of the operators act, see equations 2.10- 2.12 for these details. 

The requirement that the Gibbs-Boltzmann distribution \1/ GB given by equa- 
tion 5.1 be invariant under the stochastic dynamics is equivalent to the distribution 
yielding V • J = 0. We find it convenient to group terms and express this condition as 

(5.6) V ■ J = Ax + A 2 + V • A 3 + V • A 4 = 



16 



P.J. ATZBERGER 



where 
(5.7) 

Ai = [(A + V X A • v + Ai) • V P E + (-V x $ + Ci) ■ V v £ + (v) • V X E] (-^r) -1 ^ c 
A 2 = [V p • (A + V X A • v + Ai) + V v • (-V x $ + Ca) + V x • (v)] * GB 

A 3 = -^(V-G)* GB 



Cu + A 2 + [GppVpS + G pv V v E + G p xVx£] (2k B T)~ 1 
-T + C2 + [G vp V p E + G VV V V E + G vX V x £] (2k B T)- 1 
[G Xp V p E + GxvV v £ + GxxV x £] (2k B T)- 1 



*GB- 



We assume here that the Lagrange multipliers can be split A = Ai + A2 and C = Ci + C2 
to impose the constraints by considering in isolation different terms contributing to 
the dynamics, see equation 5.7. This is always possible for linear constraints. The 
block entries of the covariance operator G are denoted by Gij with i,j E {p, v,X}. 
For the energy of the discretized system given by equation 2.4 we have 

(5.8) V Pn £ = u(x n )Aa£ 

(5.9) V Vq E = J2 u ( x m) • (-V v ,A[mv]m) Aa&, + mv, 

m 

(5.10) V x ,£ = u ( x -) ' (-Vx,A[mv] m ) Ax 3 m + V x ,$. 

m 

where u = p _1 (p — Ajrov]). Similar expressions for the energy of the undiscretized 
formalism can be obtained by using the calculus of variations [26] . 

We now consider V • J and each term A\,Ai, A3, A4. The term A\ can be shown 
to be the time derivative of the energy A\ = dE/dt when considering only a subset 
of the contributions to the dynamics. Thus, conservation of the energy under this 
restricted dynamics would result in A\ being zero. For the SELM formalism, we 
find by direct substitution of the gradients of E given by equations 5.8- 5.10 into 
equation 5.7 that A\ = 0. When there are constraints, it is important to consider 
only admissible states (p,v,X). This shows in the inviscid and zero temperature 
limit of SELM the resulting dynamics are non-dissipative. This property imposes 
constraints on the coupling operators and can be viewed as a further motivation for 
the adjoint conditions imposed in equation 2.5. 

The term A2 gives the compressibility of the phase-space flow generated by the 
non-dissipative dynamics of the SELM formalism. The flow is generated by the vector 
field (A + VxA-v + Ai, — Vx^ + Cii v ) on the phase-space (p, v, X). When this term 
is non-zero there are important implications for the Liouville Theorem and statistical 
mechanics of the system [62]. For the current regime, we have A 2 = since in 
the divergence each component of the vector field is seen to be independent of the 
variable on which the derivative is computed. This shows in the inviscid and zero 
temperature limit of SELM, the phase-space flow is incompressible. For the reduced 
SELM descriptions, we shall see this is not always the case. 

The term A3 corresponds to fluxes arising from multiplicative features of the 
stochastic driving fields. When the covariance G has a dependence on the current state 
of the system, this can result in possible changes in the amplitude and correlations 
in the fluctuations. These changes can yield asymmetries in the stochastic dynamics 
which manifest as a net probability flux. In the SELM formalism it is found that in 
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the divergence of G each contributing entry is independent of the variable on which 
the derivative is being computed. This shows for the SELM dynamics there is no such 
probability fluxes, A3 = 0. 

The last term A4 accounts for the fluxes arising from the primarily dissipative 
dynamics and the stochastic driving fields. This term is calculated by substituting the 
gradients of the energy given by equation 5.8- 5.10 and using the choice of covariance 
structure given by equations 2.13 or 4.11. By direct substitution this term is found 
to be zero, A 4 = 0. 

This shows the invariance of the Gibbs-Boltzmann distribution under the SELM 
dynamics. This provides a rather strong validation of the stochastic driving fields 
introduced for the SELM formalism. This shows the SELM stochastic dynamics are 
consist with equilibrium statistical mechanics [54]. 

5.2. Regime II. For the reduced stochastic dynamics given by equations 2.10- 
2.12, the probability density \&(p, X, t) satisfies the Fokker-Planck equation with the 
probability flux 

*gb-^(V-G)* gb -^GW gb . 

The G denotes the covariance operator for the stochastic driving fields given by equa- 
The invariance of the Gibbs-Boltzmann distribution requires 

= A 1 + A 2 + V • A 3 + V • A 4 = 

= [(A + (V x • A)fc B T + AO • V p £ + (p-'T + Ci) • V x £] (-fc B T)- l * OB 
= [V p ■ (A + (V x • A)fc B T + Ax) + V x • (p _1 r + &)] *gb 

= -~(V-G)* OB 

" {p- l C + A 2 ) + [Gp P V p £ + GpxVx-E] (2fc B T)~ 1 1 
(T- 1 + C2) + [GxpVp£: + GxxVx£](2fc i3 T)- 1 J 

To simplify the notation we have suppressed explicitly denoting the functions on 
which the operators act, which can be inferred from equation 2.16- 2.17. In the 
current regime m = and the energy given by equation 4.10 has gradients given by 

(5.13) Vp n £ = u n Aa;3 

(5.14) Vx,-E = V x ,$. 

Similar expressions can be obtained for the undiscretized formalism using the calculus 
of variations [26] . 

We now consider V- J and Ai,A 2 , A3, A4. The terms have a similar interpretation 
as in Section 5.1. The A\ term can be interpreted as the time derivative of the energy 
A\ = dE/dt when considering only a subset of the contributions to the dynamics. 
By direct substitution of the gradients given by equation 5.13- 5.14, we find A\ = 
— ((Vx • A) • W p E)^ GB . This differs from the non-reduced equations in which this 
term was zero, see Section 5.1. 

The term A 2 gives the compressibility of the flow generated by the vector field 
(A + (Vx ■ A)fc B T + Ai,p _1 r + Ci) on the phase-space (p, X). For the reduced 
equations, the phase-space flow has compressibility given by A 2 = (p _1 Vx • r)\l' GB , 



(5.11) J = 



P - 1 C + A+ (V x • A)k B T + A 
p^T + T- 1 + C 



tion 4.12. 

(5.12) 
V J 

A 1 : 

A 2 

A 3 : 
A 4 : 



18 



P.J. ATZBERGER 



which in general is no longer zero. However, we have that A\ + A2 = 0. This follows 
from the form of the gradients given by equation 5.13- 5.14 and from the properties 
of r and A. In particular, that the operators are linear and that they are adjoints 
r = A' in the sense of equation 2.5. 

The term A3 accounts for probability fluxes driven by multiplicative features of 
the stochastic driving fields. It is found this term is zero A3 = 0. This follows from 
the divergence in which each entry of G is independent of the variable on which the 
derivative is applied. The term A4 accounts for fluxes arising from the dissipative 
contributions to the dynamics and the stochastic driving fields. By direct substitution 
of the gradients given in equation 5.13- 5.14, and the choice made for G given in 
equation 4.12, we find this term is zero, A4 = 0. This establishes V • J = and that 
the Gibbs-Boltzmann distribution is invariant for the SELM dynamics. 

5.3. Regime III. We now discuss briefly the reduced stochastic dynamics given 
by equations 2.21- 2.23, in which T — > 00. In this regime, the probability flux is 
almost identical to equation 5.11, with any terms involving T _1 set to zero. With 
this substitution, it immediately follows that the Gibbs-Boltzmann distribution is 
invariant under the SELM dynamics. 

5.4. Regime IV. In the over-damped regime in which the fluid is no longer ex- 
plicitly represented, we have the reduced stochastic dynamics given by equations 2.24- 
2.26. The probability density ^(X, t) for the current state of the system is governed 
by the Fokker-Planck equation with the probability flux 



In this regime, the covariance is given by G = 2/cbT7J selm , see Section 3.4. This gives 
Vx • G = 2k B TV x ■ if SELM and ^GV X * GB = # SE lm[-V x $]*gb. Substituting these 
expressions in equation 5.15, we find J = 0. This establishes for the over-damped 
regime of the SELM formalism that the Gibbs-Boltzmann distribution is invariant 
and satisfies detailed balance. 

5.5. Summary. For the SELM formalism, we have demonstrated in each regime 
that the Gibbs-Boltzmann distribution is invariant. This shows the SELM formalism 
yields appropriate behaviors with respect to equilibrium statistical mechanics. 

6. Applications. To demonstrate how the SELM formalism can be used in 
practice, we consider spherical particles which have translation and rotational degrees 
of freedom. We give specific operators representing the coupling of the particles and 
fluid. We compare this SELM formalism with classical results from fluid mechanics. 
We should mention that similar approaches for the SELM formalism can be applied 
much more generally to represent spatially extended structures, such as filaments, 
membranes, or even deformable bodies. The development of representations and 
specific coupling operators for these structures will be the focus of future work. 

6.1. Particles with Rotational and Translational Degrees of Freedom. 

To describe particles which can exhibit translational and rotational motions, we use 
the degrees of freedom X cm for the center of mass and for the rotational configura- 
tion. To describe the full configuration of a particle, we define the composite vector 
X = (X cm , 0). To investigate how the coupling operators capture the hydrodynamics 
of the system, it is convenient to characterize the system in Regimes III and IV of 



(5.15) 

J = [ #s EL m[-V x *] + fe B T(Vx • fl" SELM ) ] * 



GB 



^(Vx-G)* 



GB 
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Sections 2.3 and 2.4. This highlights central features of the coupling operators also 
relevant in the other regimes. Given the specific degrees of freedom of the particles, 
it is convenient to express equation 2.22 as 

(6.1) ^ = r 0U 

(6.2) f = r lU . 

The u = /o _1 p. To represent the kinematics of such particles for a given state of the 
flow field of the fluid we use 

(6.3) r u = £ ( Vo (y m - (X cm + z)) u m >£ W=Jl Aa& 

III 

(6-4) r lU =^( m (y m - (X cm + z)) (z x u m ) ) SM=R Ax 3 m . 

m 

The angle brackets denote an average over the surface of the sphere which is given by 
the quadrature 

The Wk denote the quadrature weights and the z& denote the quadrature nodes. In 
practice, we use the Lebedev quadratures [39]. 

These SELM kinematics are in fact closely related to the exact kinematics of a 
passive spherical particle expressed using the Faxen Theorem of fluid dynamics [8]. 
The expression from the Faxen Theorem corresponds to the continuum limit of the 
above expressions and when the kernel functions are replaced by Dirac (5-functions [8, 
42, 32]. 

For particles which actively exert forces on the fluid, we develop a coupling oper- 
ator for the force by using the adjoint condition given by equation 2.5. In fact, this 
condition can be interpreted as requiring the fluid-structure coupling conserve the 
energy of the system in Regime III in the inviscid and zero temperature limit. This 
is seen formally by letting C — > and T — > in equation 2.21 and computing what is 
required for dE/dt = 0. Using this condition, we obtain the fluid coupling operator by 
considering directly the adjoint condition A = I^. For the specific coupling operators 
considered for the particles, the adjoint condition of the discretized system gives 

(6.6) A (x m ) = (( 77o(x m - (X cm + z)) } S - )W=H ) F 

(6.7) Ai(x m ) = (( zr?i(x m - (X cm + z)) > S)W=fl ) x T. 

The F = — 9$/t?X cm is the total force acting on the particle and T = — <9$/<90 is 
the total torque acting on the particle. The adjoint condition holds exactly for the 
discretized system provided the same quadrature is used for both T and A. 

In the SELM formalism, the coupling operators L and A encapsulate the effective 
hydrodynamic coupling between the particles and fluid. To characterize how the pre- 
sented SELM operators represent such hydrodynamic coupling in practice, we consider 
the interactions between two spherical particles. The coupling of two spherical par- 
ticles has three modes of coupling: (i) translation-translation, (ii) rotation-rotation, 
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(a) 



(b) 



(c) 



<« ► ^ 

(d) 



Fig. 6.1. Translational and Rotational Motions of Spherical Particles. On the left is shown 
the Lebedev quadrature nodes for N = 110 which are used for computing averages on the surface 
of a sphere used in the SELM coupling operators. On the right is shown the three modalities of 
coupling for spherical particles. These are Translation- Translation (b), Rotation- Rotation (c), and 
Rotation- Translation (d). 



and (iii) rotation-translation, see Figure 6.1. To characterize these coupling modes, 
we consider the effective hydrodynamic coupling tensor (mobility) which relates ap- 
plied forces and torques to consequent motions of the particles. An effective coupling 
tensor can be determined rather naturally for the SELM formalism by considering 
Regime IV at zero temperature. In this regime, the effective coupling tensor arises in 
the dynamics 



(6.8) 
(6.9) 



d 






dt 














^RR 




- x A fe . 



F 
T 



The tilde is used throughout to distinguish the tensor components associated with 
the SELM formalism from those arising in other theories from fluid mechanics. The 
subscripts T, R indicates components related respectively to the translational and 
rotational degrees of freedom. In practice, the components of these tensors can be 
easily computed from an implementation of SELM by applying forces or torques which 
are set to on only one particle and considering the i th component of the velocity 
or angular velocity of the other particle. 

To compare the results of the proposed SELM coupling operators with classical 
results of fluid mechanics, we consider the following Rotne-Prager-Yamakawa tensors 
(RPY tensors) [55, 66, 53] 



(6.10) H TT (r) 



^(l+^)l+^(l-H^)ft T , forr>2a 
S-Kfir V 3H J 8-K^r \ r 2 J 1 — 



3r \ 
4a I 



-rr 



for r < 2a 
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(Translation-Translation) (Rotation-Rotation) (Rotation-Translation) 

1 1 1 1 1 1 1 1 1 1 0.12i 1 1 1 1 1 1 1 . 1 1 1 1 o.Z5i 1 1 1 1 1 1 1 1 




r^ r/a ; r/a 



Fig. 6.2. Effective Hydrodynamic Coupling Tensors of the SELM Operators. The components 
of the effective hydrodynamic coupling tensor associated with the SELM operators are compared with 
the Rotne-Prager-Yamakawa Tensors of fluid mechanics. All components are scaled by the reference 
mobility mt = l/6ir[iat for translation and m r = l/8np,af for rotation. The kernel functions r)o 
and rji were chosen to be the radial symmetric function shown as the inset on the left. The other 
insets show the near-field hydrodynamic interactions. 



1 



(6.11) H RR (r) = - 3(X-3ff 

(6.12) H TR (r) = — x . 

In this notation, f = r/r and we denote by fx the matrix which represents the action 
of the cross-product of a vector with f . The a denotes the effective hydrodynamic 
radius of the particle. The RPY tensors capture far-field interactions and are expected 
to be accurate physically only when the particles are sufficiently separated. They are 
not designed to capture lubrication interactions or other near-field effects [11, 1, 55, 
66]. 

For the proposed coupling operators, very good agreement is found with the RPY 
tensors for all of the different coupling modes, see Figure 6.2. This agreement is espe- 
cially good for particles separated at least a distance of two radii. At closer distances 
the RPY tensors either reflect some type of regularization or they diverge. In the 
case of translation-translation coupling the SELM formalism agrees even in the near- 
field with the RPY tensor [55, 66]. For the other cases, the RPY tensors diverge. In 
all cases, the SELM formalism provides as the separation distance becomes small a 
regularized tensor. The regularization can be interpreted as an interpolation between 
the two particle far-field interaction to the single particle response, see insets in Fig- 
ure 6.2. For the rotation-translation coupling model, we see the operators have the 
important property for the single particle response that the rotation and translational 
motions are decoupled. This agrees with predictions from fluid mechanics for linear 
Stokes flow. 

These results for the proposed operators demonstrate that the SELM formalism 
provides a practical approach for simulating hydrodynamically coupled spherical par- 
ticles. We remark that the near-field artifacts are a by-product of the specific coupling 
operators utilized to approximate the fluid-structure interactions. For many practical 
systems, repulsive long-range interactions keep particles well-separated avoiding these 
near-field artifacts. It should be emphasized the SELM formalism is not limited to 
such cases, since a more accurate approach which captures near-field effects can be 
developed by using a different choice for the coupling operators. 

While we have presented only one rather special application, the SELM formalism 
can be applied more broadly. This simply requires appropriate representations for the 
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structures and a choice for the coupling operators. It is expected a wide variety of 
structures could be studied using this approach, including particles of non-spherical 
shape, filaments, membranes, and dcformable bodies. The development of coupling 
operators for these systems will be the focus of future work. 

7. Conclusions. An approach for fluid-structure interactions subject to ther- 
mal fluctuations was presented based on a mechanical description utilizing both Eule- 
rian and Lagrangian reference frames. General conditions were established for opera- 
tors coupling these descriptions. Reduced descriptions for the stochastic dynamics of 
the fluid-structure system were developed for several physical regimes. Analysis was 
presented for each regime establishing for the SELM stochastic dynamics that the 
Gibbs-Boltzmann distribution is invariant. The SELM formalism provides a general 
framework for the development of computational methods for applications requiring 
a consistent treatment of structure mechanics, hydrodynamic coupling, and thermal 
fluctuations. 
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